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f^*^ i Abstract. An identification of a spherically symmetric potential by its phase shifts is an im- 

, portant physical problem. Recent theoretical results assure that such a potential is uniquely 

j i ■ defined by a sufficiently large subset of its phase shifts at any one fixed energy level. However, 

Q ' two different potentials can produce almost identical phase shifts. That is, the inverse problem 

i of the identification of a potential from its phase shifts at one energy level k 2 is ill-posed, and 

the reconstruction is unstable. In this paper we introduce a quantitative measure D(k) of this 
1 instability. The diameters of minimizing sets D(k) are used to study the change in the stability 

£N| ' with the change of k, and the influence of noise on the identification. They are also used in 

the stopping criterion for the nonlinear minimization method IRRS (Iterative Random Reduced 
T— H ■ Search). IRRS combines probabilistic global and deterministic local search methods and it is 

1 used for the numerical recovery of the potential by the set of its phase shifts. The results of the 

identification for noiseless as well as noise corrupted data are presented. 

m : 
o ■ 

o ■ 

i3 ' 1. Introduction 

Identification of a material from the measurements of its bombardment by particles has long been 
!-£h ■ of interest in physics. The results of such an experiments are used to determine the phase shifts. 

^p" 1 ! Identification of a potential by its phase shifts is an important physical problem, see J|, [| |[ [ll], [TJ] . 

r*| ■ It has recently been shown that a sufficiently large infinite subset of the set of single energy phase 

shifts uniquely defines a bounded compactly supported potential , ]T^ |. However, in practice, only 
a finite (and relatively small) subset of the phase shifts can be determined from experimental data. 
It has been shown [jjj ff| [l7| [l9| that given such a finite set of shifts at an energy level k 2 one 
can find several quite different potentials having practically the same phase shifts at all I (angular 
momenta) although Ramm's uniqueness theorem |l2| guarantees that the phase shifts 8i :— 5(1, fc), 

; I ■ known at a fixed k > for all values of I £ £., such that X);e£ i^o 7 = 00 > determine uniquely a 

compactly supported spherically symmetric potential q(r), such that q(r) — when r > a > and 
Jo r 2 \q(r)\ 2 dr < oo. Thus, the inverse problem of the recovery of the potential by its phase shifts 
at one fixed energy level can be severely ill-posed, Jig ]. In this paper we describe a quantitative 
measure D(k) of such an instability, and an algorithm for its computation. D(k) is the diameter 
of the minimizing set defined in Section 3. The algorithm provides a method for its numerical 
estimate. Several numerical examples are presented in Section 4. 

Let q(x), x £ R 3 , be a real- valued potential with compact support. Let R > be a number such 
that q(x) = for \x\ > R. We also assume that q £ L 2 (B R ) , B R = {x : \x\ < R, x £ R 3 }. Let S 2 
be the unit sphere, and a £ S 2 . For a given energy k > the scattering solution ip(x, a) is defined 
as the solution of 

(1.1) Alp + k 2 i) - q(x)i> = , xeM 3 

satisfying the following asymptotic condition at infinity: 
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(1.2) i^ = ^jo + v, ^ a :=e zka - x , a e S 2 



(1.3) lim 
It can be shown, that 



dv , 

— ikv 

or 



ds = 0. 



& ikr / y \ x 

(1.4) ip(x, ot) = V'o + A(a' , a, k) \- o I - I , as r —* oo , — = a' . 

r \r J r 

The function A(a',a, k) is called the scattering amplitude, a and a' are the directions of the 
incident and scattered waves, and k 2 is the energy, see ||, jyj . 

For spherically symmetric scatterers q(x) = q(r) the scattering amplitude satisfies A(a! , a, k) = 
A(a' -a, k). The converse is established in [Q. Following the scattering amplitude for q = q(r) 
can be written as 



CXD / 



(1.5) A{a',a,k) = Y, Mk)Yi m (a')Y lm (a) 



1=0 m=-l 



where Yj m arc the spherical harmonics, normalized in L (S ), and the bar denotes the complex 
conjugate. 

The fixed-energy phase shifts — it < Si < tt (Si = 5(1, k), k > is fixed) are related to Ai(k) (see 



e.g., 1 18 ) by the formula: 

(1.6) i4,(fc) = ye <ffl sin(<*i)- 

2. Phase Shifts for Piecewise-Constant Potentials 

In general, phase shifts for a spherically symmetric potentials can be computed by a variety 
of methods, e.g. by a variable phase method described in 0]. The computation involves solving 
a nonlinear ODE for each phase shift. However, if the potential is compactly supported and 
piecewise-continuous, a much simpler method described in [fjj can be used. It is summarized 
below. 

Consider a finite set of points = ro < r\ < T2 < ■ ■ ■ < r/v = R and a piecewise-constant 
potential 

(2.1) q(r) = qi, on [ri_i,r<) for i = 1, . . . , N, and q — for r > R. 

Denote k 2 := k 2 — qi, where i — 1, . . . , N, and k is some fixed positive number. Consider the 
following problem for the radial Schrodinger equation: 

(2.2) ^ + ( k z- l l+fiy i = m , to(2I + l)!!r- , -Vi(r)=l, 
which we rewrite as: 

( 2 - 3 ) ( K i- ± ^ 2 W=° 



on the interval r%—i <r <r%. On [r«— i,rj) one has the following general solution of (2.3) 
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(2.4) ipi(r) = Aijifar) + B,-nj («»»"), 

where 



(2-5) Ji(kr) = y — J l+1/2 (kr) , n z (fcr) = ^ —N l+1/2 (kr) 

and J; , TV; are the Bessel and Neumann functions. 

We assume below that does not vanish for all i. If Ki = for some i, then (2.3) has the 
solution 

(2.6) W (r) = A t r l+1 + B t r~ l , 

and our approach is still valid with obvious changes. 

From the regularity of ipi at zero one gets B\ = 0. Denote Xi = Bi/Ai, then x\ = 0. We are 
looking for the continuously differentiable solution (pi. Thus, the following interface conditions 
hold: 

Aiji(Kin) + Binifcri) = A i+1 ji(Ki + in) + B i+1 m(K i+1 ri), 

(2.7) 

l^[ A ifi( K i r i) + B t n 'i{ K i r i)} = A i+1 j' l {K i+ - L r i ) + B,- +1 n{(K i+ ir,-). 
The Wronskian W(ji(r),rii(r)) — 1, thus 

A i+ i = n' l (Ki + ir l )[Aiji(mr i ) + B l ni(n i r l )] - ■^■ni(K i+1 ri)[A i j' l (K i ri) + B l n[{K i r i )], 

(2.8) 

B l+ i = 7 r^ji(K i+1 r i )[A l j' l (K l r i ) + B^k^)] - jj( K *+l r »)[ j ' 4 -iii( /s i r i) + B l ni(K i r i )\. 
Therefore 

(2.9) (B +1 )= — {f f){B 

\B l+ iJ k 1+ i \a 21 a 22 J \Bi 

where the entries of the matrix a 1 can be written explicitly: 



(2-10) 



a\ 2 = K i+ in;(Kjri)nJ(Ki + iri) - K l n' l (K i r. l )ni(n. l+1 r i ), 
"21 = ^ij'l(^in)jl(Ki+in) - K*+ljj( K i r t)#( K »+l r «)) 
a 4 22 = Ki«K K » r 0jK K i+l r i) - K l+1 ni{Kiri)j[(K i+1 ri). 



Thus 

(2.11) a; 4+1 = 

The phase shift 5(k,l) is defined by 



Qgl + "22^ , ._ A 



(2.12) W (r) ~ sin(fcr - ^ + <5(fc, Z)) r oc , 

where Z) is the Jost function. For r > R 



(2-13) 



<M r ) = AN+xj^kr) + B N+1 ni(kr). 
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From ( [2.13 ) and the asymptotics ji(kr) ~ sin(/cr — Ztt/2), n/(fer) ~ — cos(£;r — Itt/2), r — > oo, 
one gets: 

(2.14) tan£(M) = --p^ = -zjv+i . 

Finally, the phase shifts of the potential q(r) are calculated by the formula: 

(2.15) S(k, I) = — arctanirjv+i. 

Let qo(r) be a spherically symmetric piecewise-constant potential. Fix an energy level k and 
a sufficiently large N. Let {<5(fc,0}^Li be the set of its phase shifts. Let q(r) be another such 
potential, and let {S(k, l)}[Li be the set of its phase shifts. 

The best fit to data function <I>(q, k) is defined by 

(2.16) #( g , fc)= ^'y f j-^ f )' 2 , 

The phase shifts are known to decay rapidly with I, see Ew. Thus, for sufficiently large N, the 



function $ is practically the same as the one which would use all the shifts in (2.16). The inverse 
problem of the reconstruction of the potential from its fixed-energy phase shifts is reduced to 
the minimization of the objective function <E> over an appropriate admissible set. A minimization 
algorithm for this nonlinear problem providing a stability estimate for the identification of the 
original potential qq is given in the next Section. 

3. Global and Local Minimization Methods 

We seek the potentials q(r) in the class of piecewise-constant, spherically symmetric real-valued 
functions. Let the admissible set be 

(3.1) A a dm C {(n,r2, . . . ,ru,qi,Q2, ■ ■ ■ ,£?m) : < n < R, qi ow < q m < qugh] , 

where the bounds qi ow and qhigh for the potentials, as well as the bound M on the expected 
number of layers are assumed to be known. 

A configuration (n, r2, . . . , rjvf, Qi, qi, • ■ ■ , 9m) corresponds to the potential 

(3.2) q(r) = q m , for r m _i < r < r m , 1 < m < M , 

where tq = and q(r) = for r > tm = R. 

Note, that the admissible configurations must also satisfy 

(3.3) n < r 2 < r 3 < • ■ ■ < r M ■ 

Given an initial configuration Qq £ A a( i m C M 2M , a local minimization method finds a local 
minimum near Qq. On the other hand, global minimization methods explore the entire admissible 
set in order to find a global minimum of the objective function. While the local minimization is 
usually deterministic, the majority of the global methods are probabilistic. As usual for inverse 
scattering problems, the best fit to data function $ has many local minima and points of nondif- 
ferentiability, see |t]] . In this situation a combination of global probabilistic and local deterministic 
methods proved to be successful. 

In JsJ such an algorithm (the Hybrid Stochastic-Deterministic Method) has been applied for the 
identification of small subsurface particles, given a set of surface measurements. The HSD method 
could be described as a variation of a genetic algorithm and a local search with reduction. In M 



STABLE IDENTIFICATION OF POTENTIALS 



two global search algorithms in combination with a special local search method were applied to 
the identification of piecewise-constant scatterers by acoustic type measurements. The Rinnooy 
Kan and Timmer's Multilevel Single-Linkage Method in a combination with a special Local Mini- 
mization Method has been applied to the identification of piecewise-constant spherically symmetric 
potentials by their phase shifts in 0. We have used the Reduced Random Search Method in 
to find different potentials with practically the same phase shifts. 

In this paper we use the Modified Reduced Random Search Method. The important modification 
consists of the consideration of minimizing sets and their diameters as quantitative measures of 
the stability of the minimization algorithm. 

In a pure Random Search method a batch H of L trial points is generated in A al j m using a 
uniformly distributed random variable. Then a local search is started from each of these L points. 
A local minimum with the smallest value of $ is declared to be the global one. 

In our case A a d m is a box in M 2M . The uniform random variable is called 2M times to produce 
a point in this box (after the appropriate rescaling in each dimension) . Finally, the obtained values 



of r, are rearranged in the ascending order to satisfy (3.3). 

Since the Random Search method is computationally extremely inefficient, it is modified to 
reduce the number of local searches. In the Reduced Sample Random Search method one 
uses only a certain fixed fraction 7 of the original batch of L points to proceed with the local 
searches. Typically, L = 5000 and 7 = 0.01. This reduced sample H re d of 7L points is chosen to 
contain the points with the smallest jL values of <E> among the original batch. The local searches 
are started from the points in this reduced sample. This way only the points that seem to be 
in a neighborhood of the global minimum are used for an expensive local minimization, and the 
computational time is not wasted on less promising candidates. 

Let H m i n be the jL points obtained as the result of the local minimizations (7L = 50 in our 
computations). Let S m ,„ be the subset of H m i n containing points {pi} with the smallest v^L 
(0 < v < 1, we used v = 0.1) values in H m i n . We call S m i n the minimizing set. The choice of v 
determines a representative sample of global minimizers. If all these minimizers are close to each 
other, then the objective function $ is not flat near the global minimum. That is, the method 
identifies the minimum consistently. To define this consistency in quantitative terms, let |.| be a 
norm in the admissible set. 

Let 

(3.4) D = diam(S m i n ) = max{||p, —Pj\\/d a v ■ Pi,Pj^S mm }, 

where d av is the average norm of the elements in H m in ■ The normalization by d av is introduced 
to provide comparable results for different potentials. 

Large D indicates that the found minimizers pi of <I> are far apart. In terms of the Inverse 
Problem, it means that the found solution is not stable. The diameter D is a measure of such 
instability. A detailed description of an iterative version of this algorithm is given at the end of 
this Section. First, we discuss local minimization methods. 

In our minimization algorithm the Reduced Sample Random Search method is coupled with a 
deterministic Local Minimization Method. Numerical experiments show that the objective function 
$ is relatively well behaved in this problem: while it contains many local minima and, at some 
points, $ is not differentiable, standard minimization methods work well here. A Newton-type 
method for the minimization of <E> is described in (l) . We have chosen to use a variation of Powell's 
minimization method which does not require the computation of the derivatives of the objective 
function. Such method needs a minimization routine for a one-dimensional minimization of <f>, 
which we do using a Bisection or a Golden Rule method. See ^ or |t| for a complete description 
of our method. 
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Now we can describe our Basic Local Minimization Method in R 2M , which is a modification 
of Powell's minimization method ||. It is assumed here that the starting position (configuration) 
Qo S A a dm is suppied by the procedure LMM (see below), and the entry to LMM is provided by 
the global minimization part (IRRS). 

Basic Local Minimization Method. 

1. Choose the set of directions Ui , i = 1,2,..., 2M to be the standard basis in K 2M 

« i = (0,0,... ) l,...,0), 

where 1 is in the i-th place. 

2. Save your starting configuration supplied by LMM as Qo • 

3. For each i = 1, . . . , 2M move from Qo along the line defined by ui and find the point of 
minimum Q*. This defines 2M temporary points of minima. 

4. Re-index the directions Ui, so that (for the new indices) ^(Qi) < ^{Q\) <i • • ■ > ^{Q\m) — 

HQo). 

5. For i = 1, . . . , 2M move from Qi_i along the direction Ui and find the point of minimum Q.;. 

6. Set v = Q 2 m - Qo- 

7. Move from Qo along the direction v and find the minimum. Call it Qo again. It replaces Qo 
from step 2. 

8. Repeat the above steps until a stopping criterion is satisfied. 

Note, that we use the temporary points of minima Q\ only to rearr ang e the initial directions Ui 
in a different order. The stopping criterion is the same as the one in |10| , Subroutine Powell]. 
Still another refinement of the local phase is necessary to produce a successful minimization. 



The admissible set A a d m , see (|3.l|)- (l3.3| ), belongs to a 2M dimensional minimization space K 2M . 
The dimension 2M of this space is chosen a priori to be larger than 2N, where N is the number 
of layers in the original potential. We have chosen M — 6 in our numerical experiments. However, 
since the sought potential may have fewer than M layers, we found that conducting searches in 
lower-dimensional subspaces of M 2M is essential for the local minimization phase. A variation of 
the following "reduction" procedure has also been found to be necessary in [|| for the search of 
small subsurface objects, and in for the identification of multilayered scatterers. 

If two adjacent layers in a potential have values v^i and vi and the objective function $ is 
not changed much when both layers are assigned the same value Vi (or fi-i), then these two 
layers can be replaced with just one layer occupying their place. The change in $ is controlled 
by the parameter e r . We used e r = 0.1. This value, found from numerical experiments, seems to 
provide the most consistent identification. The minimization problem becomes constrained to a 
lower dimensional subspace of K 2M and the local minimization is done in this subspacc. 



Reduction Procedure. Let e r be a positive number. 

1. Save your starting configuration Qo = (ri, r%i ■ ■ ■ , tm, t>i, ^2) • • • ; %) G A a d m and the value 
$(Qo). Let the (M + l)-st layer be L M +i = {tm < \x\ < R} and v M +i = 0. 

2. Let 2 < i < M + 1. Replace v^i in the layer Li_i by Vi. This defines a new configuration 
Qf, where the layers Li^i and Li are replaced with one new layer. Here d stands for the 
downward adjustment. Compute &(Qf) and the difference cf = ^(Qo) — Q(Qf)\. Repeat 
for each layer in the original configuration Qo. 

3. Let 1 < i < M. Replace Vi + \ in the layer Li + \ by Vi. This defines a new configuration Q", 
where the layers Li and L^+i are replaced with one new layer. Here u stands for the upward 
adjustment. Compute $(Q") and the difference c" = ^(Qo) — $(Q")|. Repeat for each layer 
in the original configuration Qo- 
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4. Find the smallest among the numbers cf and c". If this number is less than e r $(Qo), then 
implement the adjustment that produced this number. The resulting new configuration has 
one less layer than the original configuration Qq. 

5. Repeat the above steps until no further reduction in the number of layers is occurring. 
Note, that an application of the Reduction Procedure may or may not result in the actual 

reduction of the number of layers. 

Finally, the entire Local Minimization Method (LMM) consists of the following: 

Local Minimization Method (LMM). 

1. Let your starting configuration supplied by IRRS be Qq = (n, r2, . . . , tm, i>i, ■ ■ ■ , Vm) <E 

Aadm • 

2. Apply the Reduction Procedure to Qq, and obtain a reduced configuration Qq containing 
M r layers. 

3. Apply the Basic Minimization Method in A a d m (~]H 2M with the starting configuration Qq, 
and obtain a configuration Q\. 

4. Apply the Reduction Procedure to Q\, and obtain a final reduced configuration Q[. 

As we have already mentioned, LMM is used as the local phase of the global minimization. The 
global part is described as follows: 

Iterative Reduced Random Search (IRRS). (at the j—th iteration). 
Fix < 7, u,/3 < 1, e > and j max - 

1. Generate another batch H 3 of L trial points (configurations) in A a d m using a uniform random 
distribution. 

2. Reduce H 3 to the reduced sample H 3 red of 7L points by selecting the points in H 3 with the 
smallest 7X values of 

3. Apply the Local Minimization Method (LMM) starting it at each point in H 3 red , and obtain 
the set H 3 min consisting of the 7L minimizers. 

4. Combine H 3 min with obtained at the previous iteration. Let S 3 min be the set of points 
from H 3 min U H 3 min with the smallest vjL values of (Use for j = 1). 

5. Compute the diameter D 3 of S J min by D 3 = max{||pi -pk\\/d a v ■ PhPk E S min } , where d av 
is the average norm of the jL elements of H 3 nin U H^ in with the smallest values of 

6. (Stopping criterion). 

Let p E S 3 min be the point with the smallest value of <& in S 3 min (the global minimizer). 
If D 3 < e, then stop. The global minimum is p. The minimization is stable. 
If D 3 > e and D 3 < f3D 3 ' , then return to step 1, and do another iteration. (Stop, if the 
maximum number of iterations j ma x is exceeded). 

If D 3 > e and D 3 > /3D J_1 , then stop, p is the global minimum. The minimization is 
unstable. The diameter D 3 is the measure of the instability of the minimization. 
We used (3 = 0.95, e = 0.01 and j ma x = 6. The choice of these and other parameters (L — 
5000, 7 = 0.01, v — 0.1 e r — 0.1) is dictated by their meaning in the algorithm and the comparative 
performance of the program at their different values. As usual, some adjustment of parameters, 
stopping criteria etc. is needed to achieve an optimal performance of the algorithm. 

4. Numerical Results 

We studied the performance of the algorithm for 4 different potentials qi(r),i = 1, 2, 3, 4. In each 
case the following values of the parameters have been used. The radius R of the support of each qi 
was chosen to be R = 3.0. The admissible set A a d m ( |3.l[ ) was defined with M = 8. The Reduced 
Random Search parameters: L — 5000 , 7 = 0.01 , v = 0.1 , e = 0.01 , (3 = 0.95 , j ma x = 6. The 
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Table 1. Phase shifts of q%(r) for k — 9. 



1 




S(k,l) 


; 


5(k I) 


I 


5{k,l) 





-0 


95151654D- 


01 


12 


0.43249567D-01 


24 


-0.58868576D- 


03 


1 


-0 


59487863D- 


01 


13 


0.78575610D-01 


25 


-0.15074621D- 


03 


2 


-0 


30344479D- 


01 


14 


0.27082102D-01 


26 


-0.34641742D- 


04 


3 


-0 


36224576D- 


01 


15 


-0.10281029D+00 


27 


-0.71940777D- 


05 


4 





14419664D- 


01 


16 


-0.18261448D+00 


28 


-0.13582152D- 


05 


■5 


-0 


35167060D- 


01 


17 


-0.17579851D+00 


29 


-0.23433795D- 


06 


6 





38359584D- 


02 


18 


-0.12758628D+00 


30 


-0.37119174D- 


07 


7 





40280065D- 


01 


19 


-0.76312741D-01 


31 


-0.54203180D- 


08 


8 





45775379D- 


01 


20 


-0.38650348D-01 


32 


-0.73237693D- 


09 


9 





99311592D- 


01 


21 


-0.16752224D-01 








10 





93668476D- 


01 


22 


-0.62688318D-02 








11 





32078999D- 


01 


23 


-0.20460976D-02 









value e r = 0.1 was used in the Reduction Procedure (see Section 3) during the local minimization 
phase. The initial configurations were generated using a random number generator with seeds 
determined by the system time. The run time was between 30 minutes to 2 ho urs on a 333 MHz 
PC, depending on the wave number k. The number N of the shifts used in ( [2.16 ) for the formation 
of the objective function $(g) was determined from the condition S(k,Q) > 10~ 7 S(k,l). So, it was 
different for different potentials qi and different wave numbers k. The upper and lower bounds for 
the potentials qi ow and qmgh used in the definition of the admissible set A a( i m were chosen to reflect 
the presumed a priori information about the potentials. While one may attempt to deduce the 
bounds from the set of given data, i.e. the phase shifts, it turns out, that an extensive admissible 
set may lead to a nonuniqueness in the sought potential. This issue will be studied elsewhere. 
Let qi (r) be the following potential 



This potential is shown in Figure 1. The upper and lower bounds for the potential qi ow = — 5.0 
and quigh — 5.0 were used for all the wave numbers k — 3, 4, 5, 6, 7, 8, 9. 

The phase shifts 6(k, I) for k = 9 (computed as in Section 2) are shown in Table 1. 

The identification was attempted with 3 different noise levels h. The levels are h — 0.0 (no 
noise), h = 0.0001 and h — 0.001. More precisely, the noisy phase shifts Sh(k,l) were obtained 
from the exact phase shifts 6(k,l) by the formula 



'4.0 


< r < 0.3 


1.0 


0.3 < r < 1.0 


-2.0 


1.0 < r < 1.9 


3.5 


1.9 < r < 2.2 


1.0 


2.2 < r < 2.4 


0.0 


r > 2.4 



8 h (k, I) = S(k, I) + (0.5 - z) ■ h ■ S max , 

where 5 max = max{\Sh(k, l)\ : 1 = 0,1,..., N}, and z is the uniformly distributed on [0, 1] random 
variable. 

The distance dij) 1 {r),p 2 (r)) for any two potentials in step 5 of the IRRS algorithm was computed 

as 
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Table 2 . Diameters D of minimizing sets for qi (r) at different noise levels h. 



k 


h = 0.000 


h = 0.0001 


h = 0.001 


3 


0.886140 


1.062616 


0.870692 


4 


0.653900 


0.629565 


1.017525 


5 


0.636675 


0.661651 


0.754354 


6 


0.456157 


0.326852 


0.585782 


7 


0.008116 


0.011184 


0.066413 


8 


0.014142 


0.004978 


0.010448 


9 


0.007881 


0.011849 


0.022112 



Table 3. Diameters D of minimizing sets for q2(r) at different noise levels h. 



k 


h = 0.000 


h = 0.0001 


h = 0.001 


3 


1.184754 


1.865304 


1.094802 


4 


1.435954 


0.843646 


0.887367 


5 


0.700282 


0.684679 


1.032036 


6 


0.000568 


0.000990 


0.693338 


7 


0.000373 


0.001486 


0.420729 


8 


0.000294 


0.001435 


0.003349 


9 


0.000190 


0.001218 


0.005308 



d( P 1 (r),p 2 (r)) = \\p 1 (r)-p 2 (r) 



where the norm is the L2-norm in R 3 . 

The results of the identification algorithm (the diameters of the minimizing sets as the function 
of the wave number k) for the potential qi(r) are shown in Table 2 as well as in Figure 2. 

The diameter D < 0.01 indicates that the potentials in the minimizing set are, practically, 
undistinguishablc. That is, the identification is stable for k > 7 and a low noise level. 

The second potential ^(r) is defined by 



1i (r) 



2.0 < r < 0.5 

1.0 0.5 < r < 1.0 

2.0 1.0<r<1.5 

1.0 1.5<r<2.0 

0.0 r > 2.0 



This positive 4-layer potential is shown in Figure 3. The upper and lower bounds for the 
admissible set were the same as for q\ . The results of the identification are shown in Table 3 and 
in Figure 4. The identification is stable for k > 6 for low noise levels. It is becoming stable for the 
noise level h = 0.001 at k = 8. 

We have also considered the identification for small potentials. In this case the potential has to 
be identified from relatively small phase shifts. We defined <73(r) = 0.1<?2(r), and 54 (r) = O.Olq^M- 
It turns out, that the identification is comparable in quality to the ones above, provided that 
reasonable a priori bounds for the potentials are supplied. We used qi ow = —0.5 and qhigh = 0.5 
for the identification of potential 33 (r), and qi ow — —0.05 and qhigh = 0.05 for potential qi{r). The 
results of the identification are shown in Tables 4 and 5, as well as in Figures 5 and 6. 
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Table 4. Diameters D of minimizing sets for qs — 0.1q2 at different noise levels h. 



k 


h = 0.000 


h = 0.0001 


h = 0.001 


3 


0.914191 


0.958988 


1.271533 


4 


0.649307 


1.432330 


0.571411 


5 


0.258550 


0.425632 


0.599525 


6 


0.000541 


0.001361 


0.318754 


7 


0.000373 


0.211909 


0.004215 


8 


0.000295 


0.002696 


0.006563 


9 


0.000170 


0.003122 


0.007665 



Table 5. Diameters D of minimizing sets for = 0.01^2 at different noise levels h. 



k 


h = 0.000 


h = 0.0001 


h = 0.001 


3 


1.233750 


0.882006 


1.681609 


4 


1.249050 


0.891463 


0.683901 


5 


0.676837 


1.082434 


0.309618 


6 


0.000565 


0.001337 


0.136704 


7 


0.000270 


0.001288 


0.004241 


8 


0.000523 


0.001006 


0.002867 


9 


0.000358 


0.000991 


0.005288 



5. Conclusions 

Recovery of a spherically symmetric potential from its fixed-energy phase shifts is a classical 
physical problem. Recent theoretical results |L2j assure that such a potential is uniquely defined by 
a sufficiently large subset of its phase shifts at any one fixed energy level. However, two different 
potentials can produce almost identical fixed-energy phase shifts (|Q, p7|). That is, the inverse 
problem of the identification of the potential by its fixed-energy phase shifts is very unstable. In 
this paper we investigate the instability of the inversion by an Iterative Reduced Random Search 
Method (IRRS). The diameter of the minimizing set D is introduced to provide a quantitative 
measure for the instability. It also serves as the stopping criterion in the IRRS algorithm. The 
results show, that for several types of piecewise-constant spherically symmetric potentials the 
identification is becoming stable (D < 0.01) for phase shifts measured at a higher energy level. It 
is also shown that the introduction of a low noise level into the data does not significantly degrades 
the identification. This method can serve as a tool for experimentalists to determine if a particular 
set of phase shifts would produce a stable identification of the underlying potential, or a higher 
energy level should be used. 
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Figure 1. Potential qi(r). 




Figure 2. Potential ^(r - ). 




Figure 4. Diameters D of minimizing sets for qi (r) at different noise levels ft. 




Figure 6. Diameters D of minimizing sets for = 0.01(72 at different noise levels h. 



